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Abstract 

The mean field approach to the Nambu-Jona-Lasinio model is de- 
veloped systematically. Approximate mean field action is obtained, 
based on the study of divergencies in the mean field action. A spe- 
cial scalar case of the approximate motion equations is studied and 
inhomogenouous solutions are discussed. It is shown that the model 
can have inhomogeneous vacuum configurations which leads to bound 
fermionic states. 

1 INTRODUCTION 

The Nambu-Jona-Lasinio model^ has been proposed in 1961, long before 
QCD arose. It is a common belief that the model can be derived from QCD 
in a low-energy limit. An example of such derivation is presented in [2]. 
However, it relies on unproven (though plausible) assumptions about the 
structure of gluon propagator in the low-energy limit. 

In a previous paper [3] it has been shown that NJL can be derived from 
QCD using the mean field approach developed by Kondo [HIS]. It's notable 
that in this approach NJL is obtained naturally together with a cutoff pa- 
rameter A which is just the effective gluon mass generated by a dynamic 
mechanism suggested by Kondo in the works mentioned. This solves the 
regularization problem in a way that is in perfect agreement with the belief 
that non-renormalizable models like NJL can only be effective field theories. 

It has been demonstrated by Nambu and Jona-Lasinio that a simple 
model with only scalar and pseudoscalar quartic terms in its lagrangian has 
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nontrivial vacuum state with fermionic condensate. Thus, initially massless 
fermions gain mass dynamically. This result was obtained in the mean field 
approach assuming spacial homogeneity, i.e. that the mean field has a con- 
stant value over the spacetime. However it has been pointed out recently that 
there can exist inhomogeneous vacuum configurations in this model, basar, 
Dunne and Thies have done this for NJL model in 1+1 dimensions [B]. 

However, if it comes to vacuum inhomogeneities, mean field approach 
becomes tricky since it's necessary to solve the non-simplified gap equation. 
Therefore, it becomes interesting to develop some approximations that would 
be less tough but allowing to grasp vacuum inhomogeneities. This is done in 
this paper. 

In the first part of this work NJL derivation from QCD is sketched briefly. 
Then the mean field approach to the model obtained is developed systemat- 
ically. An approximate mean field action is obtained, based on the study of 
divergencies in the mean field action. Afterwards a special scalar case of the 
approximate motion equations is studied and inhomogenouous solutions are 
discussed. 

2 NJL MODEL IN QUANTUM CHROMODYNAMICS 

We start with a plain (i.e.) fiavorless SU(N) model 



The cubic gluon self-interaction term can be removed by introducing an 
antisymmetric auxiliary field S^"^ in such a way that leads to the following 
lagrangian: 



where a is an arbitrary numeric parameter and fermionic terms have been 
omitted for a moment. 





d^A^ - d^A^- ig [A^,, A^] 



(1) 




(2) 
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It's easy to demonstrate that 



+ 



J VBexp \ -i I d^xTr 



(d,A, - d,A,f + 



[A^, A,f - a-'B^^ {d,A, - d,A,) - lagB^^ [A^, A,\ - ^B' 



exp {—i I d'^x 



VBexp {—if d'^x 



— TrF^ 
2 



- -Tr {B^, + a-' {d^A, - d^A,) + iag [A^, A,] }' 



--TrF"^- 
2 



= Cexp^-iJ d^x -^TrF^ 
In order to get rid of of the fourth order gluon term we note that 



(3) 



Tr [A^, Ajf = TtA^ K, [A^ A^]] = TrA, {{A,, {A^ A''}} - {A^, {A,, A^}}) = 



TT{A,,A,Y-Tr{A,,A^y 



— Tr 



{A,. A^} - ^ {A„ A^} - - -Tr {A„ A^} 



2 3, 



Then we introduce auxihary fields Qf^^, — O^^, + r^^T'"(©'^ a — 0) 
$ = + (^"T" in order to arrive at 



C^Tr 



l + a- 



-$ {A„ A^} - <j-^B^' {d,A, - d,A,) - lagB^'' [A„ A,] - 
1 2 



--B'-- 

2 2^2(1 + ^2 



3^2(1 + ^2 



(5) 



It's obvious now that the Faddeev- Popov lagrangian for gluon fields can 
be transformed into 



(4) 



C = A^-ICliA'"' + J^{B)A'^^^ + Up + Cmf{^, e,B) + ij [iD - mo j ^ (6) 

, where Cfp is the ghost part of the lagrangian and >Cmf(^,©,-B) depends 
on auxiliary fields only. The /C operator is 
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- ^(rf"V+^"V) + f r'"^5^., (7) 

where ^ is a gauge-fixing parameter. 

The lagrangian now is quadratic in gluon fields. Thus the gluons can be 
integrated out and the effective action becomes 



Seff — —iTr In i)C + 
1 



/ 



+ / d^xTr 



1 



-6^ + 



The current is 



(8) 



(9) 



= j; {B, C) + g^^.T^^. 

where (S, C) is made of the B field and Faddeev-Popov ghosts . 

It's easy to find the gluon condensate now. Let us put -B = O,0 = O,(/7" = 
0, = 00 7^ which is the case with highest unbroken symmetry. Let us also 
fix the gauge as follows: 

1 + + r' = 0. 



Then the equation for (po follows immediately: 



(1 + a^) 



= 



Putting 



a — ip, p > 1 



(10) 

(11) 
(12) 



we obtain a nontrivial solution 



00 = ± 



Wick rotation was made here so 



(13) 
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This integral is divergent and needs regularization but we won't go into it. 

Having this solution, we can build the mean field approximation which is 
equivalent to leaving the zero order term only in the saddle-point expansion 
of the path integral that corresponds to lagrangian ([6]). The mean field 
lagrangian is 



£ = -TrN' {j]^, (p2 _ 1) ^2 + (1 - p2 + ^-1) Q^Q^ + ^^^0^] A- + 

+ i)[iD- mo) ^ + Cfp (15) 
It's clear now that our solution describes dynamic gluon mass generation: 

2 ^0 

"^A = 7 16 

— 1 

Note that the effective mass doesn't depend on the free parameter p. 

Let us fix the gauge (fTOl) in the action (fTSl) and integrate out the gluons 
omitting the Faddev-Popov ghosts. The resulting action will be 



= y d'^xifj (^id — mo j ■0 + 

d^xrf^2y^(x)7^rXx)(x| {d^ + m\y' \y)i;{y)YT''^{y) 



1 



2 (p2 - 1) 



(17) 



The Green's function in the second term of the action can be expanded 
as follows: 



{x\ 



(a^ + "4)-'l!/> = %^i: (18) 



m^ n 

^ n=0 ^ 



Therefore, we can conclude that the higher order terms will be small and 
the expansion can converge if 

Ip^I < mA = A (19) 
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This condition defines the low energy hmit of the theory. 

The zero order term of the expansion above gives a four-fermion term: 

= (20) 

Using the technique developed in [3], this term can be simplified further 
which leads to the following fermionic lagrangian: 



ijji [id - mo j i)i - 



4 (p2 - 1) A2 



(21) 



Diagonality with respect to color indices is a notable feature of this la- 
grangian. It allows us to define colorless mean fields which will be done 
below. 



6 DYNAMICS OF THE FERMIONIC MEAN FIELDS 

Let us now study the NJL model itself. We'll work within the mean field 
approach to it developed first by Nambu and Jona-Lasinio to find the vacuum 
condensate of the model. Let us start with the following path integral that 
corresponds to the lagrangian fl?Il) : 

VipVi/j exp ^ (^^^ ~ "^0^ + a 'y^pifjipilj — ijj'y^tjjtjj'j^ijj 



Where we've put 



1 



a 



/3 



4(p2_i)A2 



2N 



(22) 

(23) 
(24) 



for the sake of brevity. 
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Now we introduce one scalar, one pseudoscalar, one vector and one pseu- 
dovector auxiliary fields ^, r], v^, to get rid of the fourth-order terms. We 
arrive at the following path integral: 



Z 



'Dijj'Dip'D^'Dri'Dv'Dw exp 



^2 _ ^ 2^2) 



d X 



[id — tuq + ^ + ?77^ + V + 'j^w ) ijj— 



(25) 



Let us now define the matrix mean field to be 



(26) 



It's obvious that the terms in (??) that are quadratic in auxiliary fields 
can be understood as some bilinear function of Q, therefore we'll denote them 
simply as $2(^)- 

The fermions can be integrated out now and we get the following effective 
action for Q: 



S^ff [Q] = -iNTr In i (^id - niQ + 0^ - ^ 
The equation of motion that can be derived from it is 

1 d<^2 (^) 



d'^x'i>2{n) 



(27) 



(28) 



where Gn{x, y) is the Green's function for operator id—niQ+Q. This equation 
can be called gap equation or Gor'kov equation, based on the analogy to 
superconductivity. By putting Q = const it can be reduced to a system 
of algebraic equations on the components of Q which was in fact done by 
Nambu and Jona-Lasinio themselves in the very first article on the subject. 

Studying the inhomogeneities, however, is much more complicated since 
one has to deal with a nonlinear integral equation fl28l) . Another difficulty 
is that deriving the Green's function in the right-hand side of it encounters 
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the problem of ultraviolet divergencies. Fortunately enough, we've got the 
cutoff, but it should be applied exphcitly which is what we're going to do 
now. 

Let us expand the logarithm in the effective action fl271) . We've got 



S^jf [Q] = -iNTr In i (id - mo) - iNTr In ( 1 + -r^ o ) - — [ d^x<^2 

^ ' \ id — mo J 4:a J 

= j d^x^2{^) - iNSpSiO) j d^xn{x) + 

+ ~^^P J ct^xd'^yS{x — y)Vt{y)S{y — x)Vl{x) + ... 

, where S{x — y) is the plain Dirac propagator. 

The first term in the expansion can be rewritten as follows: 



Spm j d^xhix) = Spj ^l^^(-) (30) 
Let us now demonstrate that each term in the series can be reduced to 

oo 

d'^pd^x Sj {x, p) 

. For this, it's enough to demonstrate that 

/ d'^ySix - y)h{y) [ d^pe-'P^y-'^E{p, z) = [ d^pe-"p'^'"'%{p, z) 
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So 



J d'yS{x - y)Q{y) J d'pe-''''-y-'^E{p, z) ^ 



q'^ — uIq 



(27r)4 



q^ — ml 



X 



X 



E 

n=0 



(y - z)"K.{y - z) 



d d 



n{z) 



e-ip(y-')E{p, z) = 



dSjd^qd^p -^qfa-^) q + mp 



q^ — ml 



X 



X 



n=0 



nl 



d d 



-n{z) 



d d 



-ip{y-z) 



dpai "dpar, 



z) 



d yd qd P iq(^-y) q + mp ip(y^,) 
(27r)4 q^-mf 



d d 



n=0 

= J d^qd^pe-'^'^''-^'^ 



n{z) 



d d ^ 



dpai dpa^ 



i(ox-vz) Q + IT^O f d^y_ i{p-q)y 



q^ — ml J (27r 



X 



E 

n=0 

/ 



n! 



d d 



n{z) 



d d ^ 



dpa, dpa„ 



d^qd^pe-'^^^-P^^^^r^^iP -<l)y 
q^ — mi ^ 



n=0 



d d 



n{z) 



d d ^ 



dpai "dpa„ 



d'^pe 



-ipix-z)P + ^o ^ i-i) 



E 



d d 



Q{z) 



Therefore, it's now obvious that if we define the operator 



-E{p,z) Q.E.D {31} 



dpai dpa„ 



'^{x,p) = 



n=0 



-^ 



d d 



dz"^ "dz""- 



Q{z) 



d d 



dpai dpa„ 



(32) 
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it's possible to express the expansion (12^ as 



Seff m 



1 



/ 



d^x(S>2{n)+iN 



I 




k 



p + mo 
— ml 



\V{x,p) ol 



k=l 







(33) 



We can easily read from this that all the integrations in (1291) reduce to 
expressions of type 



First of all, it's easy to see that integrals like this are nonzero for even 
values of m only and are finite if 



. It's obvious that differentiation of the integrand in flM|) increments 
its d{n, m) by 1 while leaving the difference n — m unchanged. Thus we 
can conclude from fl33l) that the series contains only the integrals fjM|) with 
m < n. But it's also true that a k-th order in Q term can contain only the 
integrals with n > k. 

Therefore we immediately conclude that the expansion (155]) has only finite 
number of divergent terms which are 

1. first order term; 

2. second order up to second derivatives of Q; 

3. third order up to first derivatives; 

4. fourth order with no derivatives of fi. 

Regularization by cutoff will replace the divergencies with some finite 
factors of order A" in the cutoff parameter. Therefore, the finite part can be 
treated as a small correction to the divergent one. So by omitting the finite 
terms we can build an approximation of the effective action ( 1271) that has 
remarkable features. First, it's similar in structure to the Ginzburg-Landau 
functional in superconductivity, since it contains derivatives up to second 
and nonlinearities up to fourth order; second, it requires no extra conditions 
being imposed on the model to be valid. 




(34) 



d{n, m) = 2n — m — 4 > 



(35) 
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To calculate the terms of our interest we use the Passarino-Veltman re- 
duction [7] 



I{2m + l,n)o,,..a2m+i = 



/(2m, ?T,)qj Q,2^ — CiS* {l]aia2--Va2m-ia2ni,) j ^'^P' 



2m 



^p2 _ 



Tn (36) 



where C is a symmetry factor and S {i]aia2--Va2,n^ia2m) is a symmetric tensor 
power of the Lorentz metrics. 

After rather long computation we obtain the following expression for the 
approximate " Ginzburg- Landau" action: 



+ 1 



327r2 
Z(A) - 2ml 



d^xSp \ 2mo [Z{A) - ml] (x) + ml h{x 



-f^Q{x)-f^Q{x) + mon{x)idn{x)- 
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mo 



dn{x) 
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-/JdMx)]Yd'^ix)- 



n{x)-ff,n{x)-f^n{x) 



+^^^h{x)r \dMx)] - '-j^h{x)r^{x)dn{x)+ 

L J b 



+—Sp 
24 ^ 



'y^Q{x)r^{x) ' + ^Spj,Q{x)-fMx)r^{x)Y^{x) \ ,(37) 



where it has been put 



Z{A) 



A' 



In 



A2 



+ 1 



(38) 



Note that differential and nonlinear terms are of the same logarithmic order 
in A. This implies that nonlinearities play an important role in the inhomo- 
geneous case and can't be treated as just small corrections. 

The only special case that should be discussed here is the zero initial 
mass case mo = when infrared divergencies appear. This seems an obstacle 
but let us recall that this is in fact the case studied by Nambu and Jona- 
Lasinio where the vacuum condensate appears. Thus we can conclude that 
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the infrared divergencies are an artifact of expansion around the false vacuum 
f2 = and the disease can be cured by expanding around the true vacuum 
of the theory. Putting fl = ^ = const in the gap equation fl28|) we obtain 



TV 



167r2 (p2 _ 1) 



In ( 4^ + 1 

i-e 



A2 



(39) 



From this equation we can conclude that C, = tA where r is a constant 
multiplier. Let us recall the definition of I{n,m) If we expand around 

the true vacuum Q = ±M = ±tA then I{n, m) get replaced with {I){n, m) = 
M~^""'"'^"'"^S'(n, m) where 6'(n, m) doesn't depend on M and A. This can 
make our approximation even more plausible for the zero mass case than 
for nonzero mass since the finite terms get multiplied by A^^""'"™'"'"^. The 
complete approximate action for the zero-mass case can be written out by 
making the following substitutions in fl37|) : 

A2 



n^VL]^2 (^) ^ $2 (^^^ - ; Z(A) ^ Z(A, M) 



INHOMOGENEOUS SCALAR VACUUM 



It seems hard to study the action f l37|) . However it posesses a simple scalar 
sector. By putting 

h{x) = (40) 

in the corresponding equation of motion the following equation can be ob- 
tained: 

-G{N,p)Z{K)i{x) - 2mo [Z{K) - m^] + 2 [Z(A) - Sm^] ^(x) - 
-Ui{x) + 6mo^^(x) - 2i^{x) = 0, (41) 

where 



A^ 

The most interesting case here is rriQ = when it reduces to the well- 
known 0^ equation: 

□C(x) - [2 - ^(A^, p)] Z(A, M)e(x) + 2C\x) = (42) 
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This equation can be obtained by making the expansion around the true 
vacuum value of Q discussed in the previous section and then changing the 
dependent variable back to 

h = -M + n 

It has the famous one-dimensional kink solution 

^{z) = Atanh{Kz) , (43) 

which corresponds to a domain wall between the areas with different vacua. 

The equation fl^Tl) doesn't have the kink but it has a periodic solution 
that can be expressed in terms of elliptic functions which is similar to what 
was obtained by Thies et. al. It can also have a soliton solution as a special 
case. 

5 CONCLUSIONS 

We've studied the effective NJL mean field action. We've built the approxi- 
mate action fl371) that doesn't require any extra conditions to be imposed on 
the dynamic variables and is governed by the cutoff parameter only. So we 
can conclude that our approximation is valid under the same conditions as 
the NJL model itself. 

The mean field dynamic in our approximation has a nontrivial scalar sec- 
tor that is described with a variant of the 0^ equation in the zero initial mass 
case. Thus in this case we've got some nontrivial inhomogeneous vacuum 
configurations, in particular the kink solution. 

Let us recall that the mean field approach allows us to construct an 
approximate fermionic lagrangian which in our case is 

C = ilj(^id-mo + n^ip (44) 

where we should substitute an actual solution of the gap equation for Q. 
Therefore we can see that vacuum inhomogeneities can lead to interest- 
ing fenomena like fermionic bound states. Indeed, if we put niQ = and 
Q = /ztanh(2;) which is just the kink case we get exactly the Jackiw-Rebbi 
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problemjS] which is known to have the zero energy solution. We've stud- 
ied the problem and found that it actually posesses discrete spectrum En = 
a/ n — n) but the details will be given in a subsequent paper. 

The scalar sector studied here is only a small part of the picture. The 
approximation f l37|) can lead to other interesting fenomena. 

The most thrilling question that arises here is whether there can exist sta- 
ble compact 3D inhomogeneous vacuum configurations with bound fermions 
that can contribute to quark confinement. 
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